Giant formation rates of ultracold molecules via Feshbach Optimized Photoassociation 
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Ultracold molecules offer a broad variety of applications, ranging from metrology to quantum 
computing. However, forming "real" ultracold molecules, i.e. in deeply bound levels, is a very 
difficult proposition. Here, we show how photoassociation in the vicinity of a Feshbach resonance 
enhance molecular formation rates by several orders of magnitude. We illustrate this effect in 
heteronuclear systems, and find giant rate coefficients even in deeply bound levels. We also give a 
simple analytical expression for the photoassociation rates, and discuss future applications of the 
Feshbach Optimized Photoassociation, or FOPA, technique. 

PACS numbers: 32.80.Pj 



In recent years, several techniques, ranging from Stark 
decelerators to buffer-gas cooling, have been developed to 
obtain cold molecules [l| . Such molecules are interesting 
for a range of applications in metrology, high precision 
molecular spectroscopy, or quantum computing [3| . How- 
ever, forming ultracold stable molecules in deeply bound 
levels remains a challenge: most approaches give tem- 
peratures still considered hot (roughly 100 mK - 1 K). 
To reach the ultracold regime (below 1 mK), direct laser 
cooling of molecules is usually not effective due to their 
rich and complex level structure Instead, it is possi- 
ble to create molecules starting from ultracold atoms, via 
photo-association (PA) or "magneto-association" (MA). 
While PA occurs when two colliding atoms absorb a pho- 
ton to form a molecule 0] , MA takes advantage of mag- 
netically tuned Feshbach resosances [Bj]. 

Over the last decade, PA has been widely used to study 
long range molecular interactions and to probe ultracold 
gases 0, and MA to realize molecular condensates @ 
and investigate the BEC-BCS cross-over regime Q- How- 
ever, both methods usually lead to molecules in highly ex- 
cited states. According to the Franck-Condon principle, 
electronic transitions in PA occur at large interatomic 
distances, leading to molecules in high rovibrational lev- 
els that can either decay by spontaneous emission or col- 
lisional quenching. To stabilize the molecules in their 
ground potentials, one could use two-photon schemes Q, 
or excited molecular states with long-range wells that in- 
crease the probability density at short range. This latter 
solution requires the existense of double-well molecular 
potentials [4j and cannot be easily generalized. In MA, 
molecules are produced by sweeping the magnetic field 
through a Feshbach resonance, which occur when the en- 
ergy of a colliding pair of atoms matches that of a bound 
level associated to a closed channel. The molecules pro- 
duced by MA are in the uppermost states near dissocia- 
tion [f| and thus relatively extended and fragile. 

In this article, we investigate a new PA scheme which 
uses a magnetically induced Feshbach resonance Q to 
enhance the probability density at short range. This Fes- 
hbach Optimized Photoassociation (FOPA) allows tran- 
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FIG. 1: FOPA: Colliding atoms (1) interact via open (blue) 
and closed (green) channels due to hyperfine interactions. A 
Feshbach resonance occurs when a bound level (2) (green wave 
function) coincides with the continuum state (blue wave func- 
tion). A photon (wavelength A) can associate the atoms into 
a bound level v (3) of the ground state potential (red) with 
inner and outer classical turning points 7?i n and 7? out . 



sitions even to deeply bound levels (see Fig. [I]). Fesh- 
bach resonances and PA have been proposed to associate 
atoms Q and convert an atomic into a molecular BEC 
[Toj ] . However 



as opposed to previous proposals 
FOPA takes advantage of the whole wave function in 
a full quantum coupled-channel calculation, and is thus 
more general than the Franck-Condon principle. 

Feshbach resonances are commonly found in both 
homonuclear and heteronuclear systems with hyperfine 
interactions. We focus our attention on heteronuclear 
systems for which the presence of a permanent dipole mo- 
ment allows transitions from the continuum directly to a 
rovibrational level v of the ground electronic molecular 
12| (see Fig. [T|) . The correspondin 
'v ro \o- v PA ) @, 



states 



photoassoci- 
ation rate coefficient Kp A = (v le i<Tp A ) [2|, |s| depends on 
u re i, the relative velocity of the colliding pair, and on o~ PAl 
the PA cross section. The bracket stands for an average 
over the distribution of v xe \, here a Maxwell-Boltzmann 
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distribution characterized by the temperature T [13j . At 
low laser intensity / and ultracold temperatures, where 
only the s-wave contributes significantly, the maximum 
rate coefficient (neglecting saturation) is [ll] 

^pa = ~^l(^=i|£0R)l^=o)l 2 , (i) 

where Qt = (27: jiksT '/ '/i 2 ) 3 / 2 , and D(R) is the appropri- 
ate dipole moment for the transition between the initial 
\^e.i=o) and final \<j> Vt j=i) states corresponding to the s- 
wave (I = 0) continuum wave function of the colliding 
pair and the populated bound level (v, J = 1) wave func- 
tion. Here fee, h, and c, are the Boltzman and Planck 
constants, and the speed of light in vacuum, respectively. 

We determine \^ e .i = o) by solving the Hamiltonian for 
two colliding atoms in a magnetic field [5j, [l4j : 

H=£ + V C + Y,H? ■ (2) 

Here, V c = V (R)P° + VifflP 1 is the Coulomb in- 
teraction, decomposed into singlet (Vq) and triplet (Vi) 
molecular potentials, with the associated projection op- 
erator P° and P . The internal energy of atom j, 

a u) _ _ 
^j nt ~ ~~ffi~**i ' ^ e ^i ~ Inij) ' B, consists of the hyper- 
fine and Zecman contributions, respectively. Here Sj and 
ij are the electronic and nuclear spin of atom j with hy- 
pcrfinc constant and B is the magnetic field. Since 
the nuclear gyromagnctic factor 7„ is three orders of mag- 
nitude smaller than 7 e , we neglect it in our calculations. 




FIG. 2: Kp A in cm 3 /s vs. the B-field (T = 50 fxK, I = 
1 W/cm 2 ) for various levels (v, J = 1) of the LiNa X J E+ 
potential, starting from 6 Li(/ = |, m = — h) and 23 Na(/ = 
l,m = —1). Two Feshbach resonances at 1081 and 1403 
Gauss enhance the PA rates by several orders of magnitude. 

We solve for l^e^^o) by using the Mapped Fourier Grid 



method [15(, and by expanding it onto the basis con- 
structed from the hypcrfine states of both atoms, 

N 

|*e,i=o) = X>«(i2){|/i,mi>®|/a,m2>} a , (3) 

a=l 

where fj = ij + Sj is the total spin of atom j, and rrij 
its projection on the magnetic axis. Here, ip a (R) stands 
for the radial wave function associated with channel a 
labeled by the quantum numbers /j, m^; the Hamiltonian 
@ couples channels with the same total projection M = 
mi + m 2- As an example, we consider forming LiNa in 
the ground X 1 E + electronic state from atoms initially 
in the 6 Li(J = i,m = -±) and 23 Na(/ = Tm = -1) 
states (a = 1), using the potentials of Ref. [lq. Eight 
channels with total M = — I are coupled, and we found 
two Feshbach resonances at 1081 and 1403 Gauss. Fig. [5] 
displays Kp A as a function of the B-field into different 
levels (v, J = 1) at T = 50 and 1=1 W/cm 2 . Near 
a resonance, Kp A is drastically enhanced by up to five 
orders of magnitude, even for the lowest (v < 10) levels. 
For typical densities (n.Li = nNa ~ 10 12 cm -3 ) and an 
illuminated volume V of 1 mm 3 , N v = ^LinNa^^pA = 
2 x 10 6 molecules are are formed in v = at 1403 Gauss 
(neglecting back-stimulation Q). 




FIG. 3: Probability density |* e _ i=0 (7?)| 2 vs. B. As B nears a 
resonance, |*]/ ei ; = o| 2 increases sharply (truncated above 0.01). 
Examples of |SI/ e> ! = o| 2 off and on resonance (green planes at 
1200 and 1400 Gauss, respectively) are shown in Fig. [4] 

These giant formation rates can be understood by the 
sharp increase in the amplitudes of the radial wave func- 
tions ip a in the vicinity of a Feshbach resonance. In 
Fig. El we show the total probability density |^I' Ci ; = o(-R)| 2 
as a function of B. As the magnetic field B nears either 
of the resonances at 1081 and 1403 Gauss, |^ ei 2 = o(-R)| 2 
increases by several order of magnitudes (qualitatively 
the same for all channels). Fig. |4] shows the total ini- 
tial probability density |'i , ej ; = o(i?)| 2 on and off resonance 
(B = 1400 and 1200 Gauss, respectively): the main ef- 
fect of the resonance is the appearance of a large peak at 



3 



shorter distance near 40 ao (see top panel, inset). This 
peak is roughly located at the classical outer turning 
point R ou t of the bound state associated to the closed 
channel, usually one of the uppermost bound levels. This 
is apparent in the top panel, where this peak almost co- 
incides with the outer lobe of \4>v=aa,j=x (^?) | 2 , the upper- 
most bound level of X 1 S + . We also observe that the off- 
resonance probability density is very much reduced when 
compared to on-resonance, leading to a very weak overlap 
integral in Kp A . The lower panel shows the short sepa- 
ration range, where the on resonance probability density 
is much larger than the off resonance case, leading to 
a subtantial overlap integral in Kp A with deeply bound 
levels (e.g., v = or 4). We also note the more compli- 
cated nodal structure of \ 1 $ ey i = Q(R)\ 2 , a direct result of 
the hyperfine mixing of the entrance channel a = 1 with 
all other channels (a = 2, ... ,8). 
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FIG. 4: Probability density on- (black) and off- (red) reso- 
nance. The top panel shows that |$ ei ; =0 | 2 gets a peak near 
R ~ 40 a.u. on resonance (inset). The upper bound level 
v = 44 of the singlet ground state is also depicted. The off- 
resonance density is negligible for R < 50 a.u. The bottom 
panel illustrates the inner region where the more deeply bound 
target levels v are located (e.g. v — and 4). Again, |>I> 6 ;=0 | 
is sizable on resonance and negligible off-resonance. 

Analytical results are obtained with a two coupled 
channel model of reduced /i, in which the wave function 
ipi of the continuum state associated to the open channel 
1 (with potential V\) is coupled to the wave function i\) 2 
associated to the closed channel 2 (with V 2 ) [13] 

_ % 2 d? fjjA f Vi Vi, 2 \ (ih\ = E (1>i 

2n dR 2 \i/j 2 J V ^2,1 V 2 J \^ 2 ) \M 

We assume both coupling terms V\^ 2 and V2.1 to be real, 
and hx the threshold E\ of channel 1 at E = 0. If the 



(4) 



couplings were switched off, the solution for the open 
channel I would be tpi — > TOg while the closed channel 
2 would have a bound state 02 — ► 0o with energy Eq. 
A resonance occurs when E is near the ener gy E q of i/jq. 
The analytical solutions for Eq.(f?]) are then 17 1 



01 (i?) = V'rcgOR) +tan<5 irr (i?) 
i?->oo 1 



2 ix 



MR) 



cos<5 V irh 2 k 



sm(kR + S hg + 5) , (5) 
sin S 0o (R) ■ (6) 



Here 5 D g and S are the background and resonant phase 
shifts, while k = y/2jlE/h. The asymptotic regular and 

sm(kR 



irregular solutions are re 



and 



cos(kR + Shg). Finally, the width T(E) of 
the resonance may vary slowly with E. 

Here, scanning the _B-field is equivalent to scanning E, 
since the position Eq of the bound state in channel 2 is 
shifted by the Zeeman interaction. To first order in k, the 
s-wave phase shifts are related to the scatterin g le ngth a 
by tan(5 + 5 a g) = —ka, with <5 D g = —kab g and 14 1 



a bg 1 



A 



B-Bq 



(7) 



where a Dg is the background scattering length of the pair 
of atoms (which can slowly vary with B) , Bq is the po- 
sition of resonance, and A is related to T(E) [3]. Intro- 
ducing the analytical solutions into Eq.([T]) leads to 



A PA — A off 



1 + Citan(5 + C2sin(5| 5 



(8) 



where 



/V off 



Hit- 

TP 



I e" 



■\(MD\M S )\ 2 is the off- 
= 0) with 0i, the hnal 
T ) / (ip v \D\ip Tee ) relates to 



Qt 

resonance rate coefficient (5 
(target) state, C\ = (0„|£)|0i, 
the open channel 1, while the coupling to the bound 
state ipo in the closed channel 2 is given by C 2 = 
-^/2/nT(MD\M/(^v\D\M g ). 

The relative importance of C\ and C 2 depends on the 
nodal structure of ip v , ipreg, 0irr, and ipo- Unless i? ut 
of ip v accidentally coincides with a node in if) ie g or ^i rr , 
the overlap integral of ip v with both rcg and 0j rr are of 
the same order, leading to \C\ \ ~ 1. The relative size of 
C 2 can be controlled by the target level v. For a deeply 
bound level, i?out is at short separation where the overlap 
with tpreg is small while the overlap with 0o can be sub- 
stantial leading to \C 2 \ ^> \C\\. For very extended levels 
v, i?out of 0„ is at large separation and the overlap with 
■00 less important, leading to \C 2 \ <C \C\\. Naturally, 
these behaviours might differ for specific levels v. 

The generalization of Eq.® to several coupled chan- 
nels is straightforward as we simply add their contribu- 
tions. Furthermore, we find that only two or three chan- 
nels contribute significantly to give these giant formation 
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rate coefficients. In Fig. [5j we show Kp A for the ground 
vibrational level v — with the same parameters used in 
Fig. O The top panel depicts the scattering length a with 
the two Feshbach resonances and its analytical fit. The 
bottom panel compares the exact numerical results using 
eight coupled channels with the simple expression ([5]) . In 
both cases, the agreement is impressive. We verified that 
similar agreement was obtainable for other levels v, indi- 
cating the broad and general validity of Eq. |(HJ) . 




B (Gauss) 

FIG. 5: In (a): scattering length a for full coupled problem 
(circles) and the fit using a = a hs (B)(l — B Al Bi — B A %^ ) with 

a bg (B) = ag£ + oi(S + B ) + a 2 (B + B ) 2 . In (b): A'£ A for 
v = at 50 fiK and 1 W/cm 2 (circles) and the simple formula 
(|8]l using tan 5 = ka,b s (B)( B Al Bl + "b^bj)- The numerical 
parameters are given in each plot [l8l |. 

In conclusion, we showed that it is possible to use Fesh- 
bach Optimized Photoassociation (FOPA) to form "real" 
ultracold molecules, i.e. in deeply bound levels, in large 
quantities. In fact, the rate coefficient increases by sev- 
eral orders of magnitude, leading to giant formation rates 
near Feshbach resonances. We applied this concept to 
LiNa, an heteronuclear system with a dipole moment. In 
addition, we gave a simple analytical model describing 
the FOPA technique. As opposed to other proposals us- 
ing the Franck-Condon principle for transition near the 
turning point R ou t of the closed channel, FOPA takes 
advantage of the full wave function and its amplification 
in the vicinity of a Feshbach resonance, making it a gen- 
eral technique. In fact, FOPA could be used to do the 
spectroscopy of more deeply bound levels of excited elec- 
tronic states which are usually not reachable by standard 
PA, bridging the gap between traditional spectroscopy 
for deep levels and PA of high-lying levels realized with 
ultracold atoms. Also, by targeting levels v for which 
C\ or Ci is dominant, it is possible to determine the pa- 
rameters of the scattering length (etb g , A, and Bq) by 
pure spectrocopic measurements. This offers an accurate 
non-destructive method to first detect a Feshbach reso- 



nance, and then obtain the scattering length parameters. 
Finally, we note that this enhancement will be present 
in other manifestation of Feshbach resonances, such as 
those obtained via electric fields or magnetic dipo- 
lar interactions (e.g., in Cr j20l|). This is a very general 
technique which can be applied to bosonic, fermionic or 
mixed species, where Feshbach resonances exist. 
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